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ABSTRACT 

We present new Period- Age (PA) and Period- Age- Color (PAC) relations for 
fundamental and first overtone classical Cepheids. Current predictions rely on 
homogeneous sets of evolutionary and pulsation models covering a broad range of 
stellar masses and chemical compositions. We found that PA and PAC relations 
present a mild dependence upon metal content. Moreover, the use of different 
PA and PAC relation for fundamental and first overtone Cepheids improves the 
accuracy of age estimates in the short-period (logP < 1) range (old Cepheids), 
because they present smaller intrinsic dispersions. At the same time, the use of 
the PAC relations improves the accuracy in the long-period (logP > 1) range 
(young Cepheids), since they account for the position of individual objects inside 
the instability strip. 

We performed a detailed comparison between evolutionary and pulsation ages 
for a sizable sample of LMC (15) and SMC (12) clusters which host at least two 
Cepheids. In order to avoid deceptive uncertainties in the photometric abso- 
lute zero-point, we adopted the homogeneous set of B,V,I data for clusters and 
Cepheids collected by OGLE. We also adopted the same reddening scale. The 
different age estimates agree at the level of 20% for LMC clusters and of 10% for 
SMC clusters. We also performed the same comparison for two Galactic clusters 
(NGC6067, NGC7790) and the difference in age is smaller than 20%. 

These findings support the use of PA and PAC relations to supply accurate 
estimates of individual stellar ages in the Galaxy and in external Galaxies. The 
main advantage of this approach is its independence from the distance. 
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1. Introduction 

Classical Cepheids are crucial objects in several long-standing astrophysical problems. 
These variable stars, due to their intrinsic luminosity and regular luminosity variation, can 
easily be detected in the outskirt of the Galactic disk and in external galaxies. Cepheids are 
the most popular primary distance indicators, since the pulsation periods and the colors are 
connected with the luminosity, and therefore the Period-Luminosity (PL) and the Period- 
Luminosity- Color (PLC) relation can be adopted for estimating cosmic distances (Feast 
1999; Udalski et al. 1999a,b,c; Freedman et al. 2001; Saha et al. 2001; Fouque et al. 2003; 
Tammann, Sandage, & Reindl 2003; Storm et al. 2004). 

Classical Cepheids cover the whole range of intermediate-mass stars, and thus it is 
not surprising that these objects are typically adopted to trace Population I stars in the 
Galactic disk (Kraft & Schmidt 1963), and star- forming regions in extragalactic systems 
(Elmegreen & Efremov 1996). The use of Cepheids as tracers of young stellar populations was 
supplemented by the evidence that, if these objects obey to a PL relation, they also obey to 
a Period- Age (PA) relation. Therefore, by taking into account the Mass-Luminosity relation 
and the stellar ages predicted by evolutionary models, it was soon recognized (Kippenhahn 
& Smith 1969; Meyer-Hofmeister 1969) that an increase in the pulsation period implies 
an increase in the stellar mass, and in turn a decrease in the Cepheid age. On the basis of 
these plane physical arguments Efremov (1978) derived an empirical PA relation by adopting 
Cepheids in Galactic, in M31, and in Large Magellanic Cloud (LMC) clusters, whose age 
was independently estimated (see also Tsvetkov 1989). More recently, Magnier et al. (1997, 
hereinafter M97) derived a new semi-empirical PA relation and used Cepheids in NGC206, 
the super-association in M31, to trace the age distribution, and in turn the star formation 
history in this region, located at the intersection of two spiral arms. A similar approach was 
also adopted by Efremov & Elmegreen (1998, hereinafter EE98) and by Grebel & Brandner 
(1998, hereinafter GB98), who also derived new PA relations on the basis of a larger sample 
of cluster Cepheids. By relying on these new relations, and on the large photometric database 
current available for LMC Cepheids -more than 1000 objects- they provided useful constraints 
on the star formation history of LMC over the last 200 Myr. 

The main advantages in using the PA relation for estimating the Cepheid ages are the 
following: i) the age estimates rely upon an observable -the period- which is marginally 
affected by systematic errors such as reddening, distance, and photometric calibrations, ii) 
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The PA relation can be applied to individual objects. Therefore, relative age estimates, 
based on this method, can supply accurate constraints on the occurrence of age gradients, 
since field Cepheids outnumber open clusters, Hi) The application of the PA relation to 
cluster Cepheids supply the unique opportunity to estimate the age of the parent cluster, 
even though the photometry of main sequence turn-off stars is lacking or uncertain. However, 
empirical PA relations are affected by three drawbacks: i) current PA relations have been 
calibrated by using cluster Cepheids. However, both the zero-point and the slope of these 
relations are affected by uncertainties in the distance modulus and in the reddening of selected 
clusters, as well as by the period range covered by cluster Cepheids. ii) Cepheid PA relations 
have only been derived for fundamental Cepheids. The possibility to derive independent PA 
relations for fundamental and first overtone Cepheids would improve the accuracy of the PA 
relation in the short-period (logP < 1) range, because the intrinsic scatter is smaller. Hi) 
Cepheid ages based on the PA relations rely on the statistical assumption that the Cepheid 
instability strip has a negligible temperature width. This is a plausible assumption in the 
short-period range, but the ages of long-period (logP > 1) Cepheids should be estimated 
on the basis of a Period- Age- Color (PAC) relation. This relation accounts for individual 
Cepheid positions inside the instability strip, and therefore it is also marginally affected 
by the spread in luminosity due to evolutionary effects (M97). The reader interested in a 
detailed discussion concerning the width in temperature of the instability strip is referred to 
Bono et al. (2000a) and to Bono, Caputo, & Marconi (2001). 

The main aim of this investigation is to supply a new theoretical calibration of both 
PA and PAC relations for fundamental and first overtone Cepheids, by adopting homoge- 
neous and detailed sets of nonlinear, convective pulsation models together with up-to-date 
evolutionary models. In order to constrain the accuracy of the new relations, we also plan 
to perform a detailed comparison between Cepheid pulsation ages and ages based on the fit 
of isochrones to turn-off stars of the parent cluster. This test will be performed on a sizable 
sample of Magellanic clusters hosting more than two Cepheids and two Galactic clusters. 

In section 2, we describe the evolutionary and pulsational frameworks adopted in this 
investigation together with the procedures for computing current PA and PAC relations. 
In section 3, we present photometric data for cluster and Cepheid stars, and discuss the 
comparison between cluster and Cepheid ages. Finally, in section 4, we summarize our 
results and briefly outline future developments. 
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2. Theoretical models 

2.1. Stellar evolutionary models 

The current evolutionary framework is based on the updated set of stellar evolutionary 
models recently computed by Pietrinferni et al. (2004). The reader interested in a detailed 
discussion concerning the input physics is referred to this paper. The main changes when 
compared with our previous investigation on intermediate-mass evolutionary models (Bono 
et al. 2000b) are the following: 

• Opacity - Conductive opacity for electron degenerate matter was estimated following 
the prescriptions by Potekhin (1999). 

• Energy loss rates - The energy loss rates for plasma-neutrino processes have been 
updated by using the most recent and accurate estimates provided by Haft, Raffelt & 
Weiss (1994). For the other processes, we follow the prescriptions adopted by Cassisi 
& Salaris (1997). 

• Nuclear reaction rates - The nuclear reaction rates have been updated by using the 
NACRE database (Angulo et al. 1999), with the exception of the 12 C(a, 7) 16 reaction. 
For this reaction we adopt the recent determination by Kunz et al. (2002). Electron 
screening is treated according to Graboske et al. (1973); 

• EOS - The updated Equation of State (EOS) by A. Irwin 1 has been used. It is 
noteworthy that this EOS allowed us to compute self-consistent stellar models in all 
evolutionary phases relevant for this investigation. For a thorough discussion concern- 
ing the key features of this EOS, see e.g. Cassisi, Salaris, & Irwin (2003) and Irwin et 
al. (2004, in preparation). 

• Convection - Superadiabatic convection is handled according to the Cox & Giuli 
(1968) formalism of the mixing length theory (Bohm-Vitense 1958). The mixing length 
parameter has been fixed according to the solar calibration (ml = 1.713). 

The efficiency of the physical mechanisms that govern the size of the convective core 
during central H-burning phases is still an open problem (see, e.g., Testa et al. 1999; 
Maeder & Meynet 2000; Barmina, Girardi, & Chiosi 2002; Brocato et al. 2003). 
Current stellar models do not account for core overshooting. The reader interested in 



1 The EOS code has been made available at ftp://astroftp.phys.uvic.ca under the GNU General Public 
License (GPL) 



- 5 - 



a detailed discussion concerning the impact of mixing processes in intermediate-mass 
stars is referred to Castellani et al. (1985), Chiosi & Maeder (1986), and Cassisi (2004). 

• Chemical composition - The sets of evolutionary models have been computed for 
Z=0.004, 0.008, and Z=0.0198 (solar). The initial He abundance was chosen according 
to Salaris et al. (2004) equal to Y = 0.245. To reproduce the calibrated initial solar 
He- abundance we adopted dY/dZ 1.4 (Pietrinferni et al. 2004). The adopted He 
abundances for the three selected metal abundances are: Y=0.251, 0.256, and 0.273. 

• Stellar mass - The mass values of the evolutionary models adopted in this investigation 
range from 2.3M Q to 1OM . The mass loss is treated according to the Reimers formula 
(Reimers 1975) with r] equal to 0.4 2 . 

Fig. 1 shows the H-R diagram for stellar models adopted in this investigation. Stellar 
isochrones have been computed with steps of 50 Myr for ages older than 200 Myr and of 20 
Myr for younger ages. Theoretical predictions have been transformed into the observational 
plane by adopting the bolometric corrections and the color-temperature transformations 
provided by Pietrinferni et al. (2004). The current set of stellar isochrones covers the 
typical age range of intermediate-mass stars with a very good resolution. Therefore, we also 
estimated the dependence of age on intrinsic parameters. In particular, we found: 

logt = (8.911 ± 0.018) + (0.423 ± 0.002) x M™ + (0.200 ± 0.008) x logZ r 2 = 0.99 [1] 

where t is the age in yr, My the TO visual magnitude, r the correlation parameter, and 
the other symbols have their usual meaning. Note that current analytical relation applies 
to stellar structures with ages ranging from 30 Myr to 1 Gyr, and for the selected chemical 
compositions (see Table 5). 

2.2. Pulsation models 

We adopted the sequences of nonlinear, convective pulsation models constructed by 
Bono et al. (1999, 2000a,2002a) to estimate the PA and the PAC relations. This set of models 
was computed by adopting chemical compositions typical of Galactic (Y=0.28, Z=0.02) and 
Magellanic Cepheids (Y=0.25, Z=0.004, Small Magellanic Cloud [SMC]; Y=0.25, Z=0.008, 



2 The entire evolutionary database can be found at the following URL address: 
http://www.te.astro.it/BASTI/index.php. 
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Large Magellanic Cloud [LMC]) and cover a wide range of stellar masses and effective tem- 
peratures. Note that the He abundances adopted to construct the pulsation models present 
a mild difference when compared with the evolutionary ones (see §2.1). However, such a 
difference has a marginal impact upon pulsation properties, and in turn on the conclu- 
sions of our current investigation. The input parameters adopted for fundamental and first 
overtone models are listed in Table 1. Both canonical and noncanonical Cepheid models 
were constructed by adopting the same mass-luminosity (M/L) relations used by Bono et 
al. (1999). The grid of canonical Cepheid models was implemented with new sequences 
for M/M = 9.4,logL/L = 4.4, Z = 0.004 and for M/M = 10.1, log L/L Q = 4.4 and 
Z = 0.008. On the basis of these models, we performed a linear least-squares fit to lumi- 
nosities and effective temperatures, i.e. logT e = a + /31ogL/L , of blue and red edges of 
both fundamental and first overtone instability strip. The analytical relations for the three 
different chemical compositions are listed in Table 2. Note that the periods of fundamental 
pulsators range from logP = 0.17 to logP = 2.02, while for first overtones they range from 
logP = —0.20 to logP = 0.88. This Table also gives the analytical relations for the effective 
temperature of the fundamental instability strip at solar chemical composition provided by 
Petroni et al. (2003). These models have been constructed by adopting the same treatment 
for turbulent convection, but with different assumptions concerning input physics and spa- 
tial resolution across the ionization regions. For these models a M/L relation based on more 
recent canonical evolutionary models was also adopted. A glance at the two different sets of 
predictions (see Table 2) indicates that the difference is at most of the order of 0.02 dex in the 
zero-point and smaller than 0.01 dex in the slope. This finding suggests that the temperature 
of the edges marginally depends on these parameters. However, recent empirical evidence 
(Tammann, Sandage, & Reindl 2003; Kanbur & Ngeow 2004; Sandage, Tammann, & Reindl 
2004) suggests that the PL relations for fundamental Cepehids, and in turn, the edges of the 
Cepheid instability strip show a change in the slope at logP « I. According to theoretical 
predictions based on convective models, this occurrence is caused by a twofold effect. In the 
long-period range (logP > 1.5) the instability strip moves toward redder colors causing a 
flattening in the PL relation (Bono & Marconi 1998; Bono et al. 2000a). This behavior is 
supported by observations (Sasselov et al. 1997). In the short-period range (logP < 0.4), 
the slope of the fundamental red edge becomes shallower, and in turn, the width in temper- 
ature of the instability strip becomes narrower (see Fig. 2 in Bono et al. 2001). This trend 
is also supported by observations (Bauer et al. 1999; Udalski et al. 1999a,b). To account 
for these effects, we also computed quadratic analytical relations connecting the effective 
temperature of the instability edges to the stellar luminosity (see Table 3). The quadratic 
fit has only been performed for canonical models, since for this set we computed a finer grid 
of models. Note that the difference between linear and quadratic relations is at most 0.05 
dex, and therefore the impact on the PA and the PAC relations is negligible. Obviously, the 
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difference increases as soon as we move from the theoretical to the observational plane, due 
to changes in the bolometric correction. We did not estimate quadratic relations for first 
overtones, because both theory and observations suggest that the instability edges for these 
pulsators are linear. The difference between fundamental and first overtone Cepheids is due 
to the narrower period and color range covered by the latter pulsators. 

The current analytical relations and evolutionary tracks (see §2.1) were adopted to select 
the evolutionary phases falling inside the instability strip. This means that we selected not 
only the second and the third crossing, but also the first one. Finally, the PA relations were 
estimated as a linear fit by weighing each individual evolutionary phase for the amount of 
time spent inside the instability strip. The coefficients of the PA relations as a function 
of the chemical composition are listed in Table 4. Fig. 2 shows the new PA relations for 
fundamental (bottom) and first overtone (top) Cepheid models. Data plotted in the bottom 
panel disclose that the dependence on the metal content is strictly decreasing when moving 
from short (old) to long-period (young) Cepheids. On the other hand, data for first overtone 
pulsators (top panel) display a different behavior for periods longer than logP > 0.5. This 
change is due to the narrowing in temperature of the first overtone instability strip when 
moving toward higher luminosities (Bono et al. 2002b). This prediction agrees quite well 
with observations, and indeed the period distribution of the first overtone pulsators in the 
Magellanic Clouds shows a well-defined cut-off toward longer periods (Udalski et al. 1999a,b; 
Beaulieu & Marquette 2000). This region of the instability strip is a robust observable, since 
it can be adopted to constrain the luminosity above which Cepheids only pulsate in the 
fundamental mode (Bono et al. 2002b). Moreover, the metallicity effect among the different 
first overtone PA relations appears negligible, within the intrinsic dispersions. 

Fig. 3 shows the comparison between current predictions for solar chemical composition 
and similar relations available in the literature. The PA relations plotted in this figure show 
a reasonable agreement, in particular if we account for the fact that they have been derived 
using different approaches to calibrate the absolute zero-point, as well as different sets of 
evolutionary models. As expected, our PA relation supplies, at fixed period, younger ages 
because the evolutionary tracks we adopted do not account for convective overshooting. 
To estimate on a quantitative basis the age difference between evolutionary tracks which 
either account for or neglect convective core overshooting, the current set of stellar models 
has been implemented with an additional set of evolutionary tracks with stellar masses 
ranging from 4 to 10 M & . These tracks were constructed by adopting a significant amount 
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of core overshooting (A c = 0.2H P 3 ) and A e (envelope overshooting 4 ) ranging from to 0.5 H p 
during hydrogen and helium burning phases (Pietrinferni et al. 2004). For each given mass 
value and chemical composition, we estimated a mean luminosity and a mean temperature 
typical of the instability strip, i.e. 3.81 < logT e < 3.85, according to the shape of the blue 
loop. On the basis of the predicted stellar mass, luminosity, and effective temperature, we 
estimated a mean fundamental period for both canonical and noncanonical models using 
the pulsation relations provided by Bono et al. (2000a). Data plotted in Fig. 4 (see also 
Table 5) indicate that the difference in age for the selected stellar masses ranges from 7% 
(8M , Z=0.008, Y=0.25) to 18% (1OM , Z=0.008, Y=0.25). Note that for M = 4M , 
Z=0.004 and M = 5M Q , Z=0.02 the noncanonical evolutionary tracks do not perform the 
loop. Interestingly enough, these are the typical differences between current PA relations and 
those based on evolutionary models which account for convective core overshooting (M97; 
GB98). 

To further improve the accuracy of Cepheid age estimates, in particular in the long- 
period range, we computed new analytical relations accounting for the width in temperature 
of the instability strip. The coefficients of the PAC relations are listed in Table 6 as a 
function of the chemical composition. Note that current PAC relations have been computed, 
for each given chemical composition, for optical colors, namely B-V and V-I. Fig. 5 shows 
the projection onto a plane of the new theoretical PAC relations. To estimate the difference 
between Cepheid ages based on the PA and on the PAC relation we computed the age 
difference for each given composition along two isoperiodic lines at logP = 1 and logP = 1.5. 
Table 7 gives for the two selected periods the age based on the PA relation (columns 2,5,8), 
and the age based on the PAC relation along the blue (columns 3,6,9) and the red (columns 
4,7,10) edge of the instability strip. Table 7 also gives the logarithmic temperatures (number 
in parentheses) of the instability edges in which the period is equal to log P = 1 or 1.5. Data 
listed in this table indicate that the age difference ranges from A logt = 0.04 — 0.07 for hotter 
Cepheids to Alogt = —0.01 0.04 for cooler ones. 



3 The extent of the non-canonical mixed region is usually denned in terms of a parameter, A, which fixes 
the length, in units of the local pressure scale height H p , covered by convective cells into the stable region 
surrounding the convective core. 

4 It is worth noting that envelope overshooting at odds with core overshooting, has no effect on the 
evolutionary lifetimes, but it substantially affects the morphology of the blue loops during the core He- 
burning phases. 
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3. Empirical data 

In order to validate current PA and PAC relations, we decided to perform a detailed 
comparison between age estimates based on the canonical Main-Sequence fitting, and on 
periods and colors of cluster Cepheids. Photometric data collected by OGLE for Magellanic 
cluster Cepheids (Udalski et al. 1999a,b; Pietrzynski & Udalski 1999) are very useful to 
accomplish this goal. These data are very homogeneous, since they have been collected with 
the same observational equipment, and reduced by adopting the same approach. Moreover, 
and even more importantly, static and variable stars have the same absolute photometric 
zero-point. To perform the test we selected LMC and SMC clusters hosting at least two 
cluster Cepheids. We ended up with a sample of 15 and 12 clusters respectively. 

Figures 6a,d and 7a,c 5 show the fit of stellar isochrones at fixed metal-content with 
cluster stars in two different Color-Magnitude Diagrams, namely V,B-V and V,V-I. Circles 
and squares mark fundamental and first overtone cluster Cepheids, according to the classi- 
fication and the membership provided by Pietrzynski & Udalski (1999). Individual cluster 
age estimates have been performed by assuming the same distance modulus (18.5 for LMC 
clusters and 19.00 for SMC clusters), and the mean reddening estimate based on cluster 
Cepheids provided by Pietrzynski & Udalski (1999). It is worth mentioning that we did not 
force stellar isochrones to fit both magnitude and colors of classical Cepheids. The extent 
in color of the so-called blue loop depends on several physical parameters and assumptions 
adopted to construct intermediate-mass evolutionary models (Chiosi, Bertelli, & Bressan 
1992; Stothers & Chin 1993,1994; Bono et al. 2000b). As a whole, the agreement between 
the fit in the V,B-V and in the V,V-I planes is quite good. The color excess in V-I was 
estimated from E(B-V) according to Cardelli et al. (1988). Table 8 lists the new cluster age 
estimates, together with the mean reddening value. To provide a plausible estimate of the 
formal error on isochrone ages, we assumed an uncertainty of 0.15 mag on the distance that 
accounts for uncertainty on the individual cluster distances and on reddening estimates and 
of 0.2 dex on the cluster metallicity. By accounting for these uncertainties and the intrinsic 
errors on the coefficients of equation 1, we end up with a tipical uncertainty on the cluster 
age of the order of 0.075 dex. 

At the same time, we estimated the cluster ages using both the PA and the PAC rela- 
tions for fundamental and first overtone Cepheids. Individual cluster ages are listed in Table 
8, together with the formal uncertainties given by the uncertainties on the coefficients of PA 
and PAC relations. Figures 8 and 9 display the relative age difference between evolution- 
ary (isochrones) and pulsation ages for LMC and SMC clusters, respectively. Interestingly 



5 Thc figures 6b,c,d and 7b, c are only available in the on-line edition of the manuscript. 



-10- 



enough, data plotted in these figures show that, within current empirical and theoretical 
uncertainties (see error bars), evolutionary and pulsation cluster ages agree on average at 
the level of 20% for LMC clusters and 10% for SMC clusters. The difference between the two 
sets, and in particular the mild increase in the LMC Cepheid ages based on the PAC relation 
(top panel Fig. 8), might be due to uncertainties in the reddening correction (Zaritsky et 
al. 2002, and references therein) and/or in the cluster membership, since LMC clusters are 
typically located across the bar. 

In order to provide a comprehensive analysis over a broad metallicity range, we also 
investigated two Galactic clusters hosting at least two classical Cepheids, namely NGC 6067 
and NGC 7790. By adopting periods, mean magnitudes, and individual reddenings provided 
by Berdnikov (2000) for Cepheids in NGC 7790 (CEa Cas, CEb Cas, CF Cas) and by Laney 
& Stobie (1994) for Cepheids in NGC 6067 (QZ Nor, V430 Nor) we found that mean cluster 
ages are logt = 7.85 ± 0.02 yr (PA) and 7.81 ± 0.02 yr (PAC) for NGC 7790, as well as 
log* = 7.76 ±0.22 yr (PA) and 7.79 ±0.11 yr (PAC) for NGC 6067 (see Table 9). It is worth 
noting that Romeo et al. (1989), by using a period-age relation based on evolutionary tracks 
accounting for convective core overshooting and mass loss, estimated an age of logt = 8.1 
yr for Cepheids in NGC 7790. 

As far as the cluster photometry is concerned, we adopted the B,V,I data available in 
the web site developed and maintained by J.-C. Mermilliod for stellar systems 
(http:/ /obswww.unige.ch/webda/webda.html). Original B,V data for NGC 6067 have been 
collected by Walker & Coulson (1985), while V,I data have been collected by Piatti, Claria, 
& Bica (1998). At the same time, B,V,I data NGC 7790, have been collected by Phelps, & 
Janes (1994), Gupta et al. (2000), and Henden (2003, 

ftp://ftp.nofs.navy.mil/pub/outgoing/aah/sequence/). Fig. 10 shows the fit between empir- 
ical data and stellar isochrones at solar chemical composition. The cluster distances adopted 
to perform the fit were estimated using the recent calibration of the optical PL relations for 
Galactic Cepheids, recently provided by Storm et al (2004), and the apparent magnitudes 
and reddenings listed in Table 9. We found that the true distance moduli are 11.27 ± 0.35 
for NGC 6067 and 12.50 ±0.23 for NGC 7790. By adopting these distances and mean cluster 
reddenings of 0.32 (NGC 6067, Coulson & Caldwell 1985) and 0.53 (NGC 7790, Lee & Lee 
1999), we estimated cluster ages of 70 ± 10 and 80 ± 10 Myr, respectively. Once again, the 
difference between pulsational and evolutionary ages is smaller than 20%. 
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4. Summary and final remarks 

We computed new PA and PAC relations for fundamental and first overtone Cepheids 
for chemical compositions typical of Magellanic and Galactic Cepheids. The current pre- 
dictions rely on homogeneous sets of stellar isochrones and pulsation models covering the 
entire range of intermediate-mass stars. We found that the metallicity has a marginal effect 
on the PA relation. To validate current theoretical scenario, we performed a detailed com- 
parison between evolutionary (isochrone fit) and pulsation (PA and PAC relations) ages of 
a sizable sample of LMC and SMC clusters hosting at least two Cepheids. In order to avoid 
deceptive uncertainties in the photometric calibration of both Cepheids and cluster stars, we 
adopted the homogeneous B,V,I photometric catalog provided by the OGLE collaboration 
(Pietrzynski & Udalski 1999). The isochrone fit was performed both in the V,B-V and in 
the V,V-I CMDs by adopting the reddening estimated by Pietrzynski & Udalski (1999) for 
cluster Cepheids. We found that pulsational and evolutionary cluster ages agree quite well, 
and indeed the difference in age ranges from 10% to 20%. We performed the same compar- 
ison for two Galactic clusters, namely NGC6067 and NGC7079. We found once again that 
the difference between pulsation and evolutionary ages is smaller than 20%. 

These findings support the evidence that PA and PAC relations supply accurate indi- 
vidual age estimates. This means that they can be adopted to trace the age distribution 
not only in the Galactic disk, but also across the main body of the Magellanic Clouds. It is 
worth mentioning that the use of independent relations for fundamental and first overtone 
Cepheids will supply more accurate age estimates in the short-period (logP < 1) range, 
since the intrinsic dispersion is smaller. At the same time, the use of the PAC relations 
will supply more accurate age estimates in the long-period (logP > 1) range. Note that 
the vertical error bars in Fig. 8 and 9 mainly depend on the number of cluster Cepheids. 
Individual Cepheid ages based on the PAC relation are more accurate than the ages based on 
the PA relation, as soon as the uncertainty on the reddening correction is on average smaller 
than 0.08 mag for B-V colors and 0.07 mag for V-I colors. These limits have been estimated 
by accounting for the color coefficients listed in Table 6, and the age differences between 
PA and PAC relations listed in Table 7. The theoretical framework we developed can be 
adopted to detect the occurrence of age gradient (s), not only across star-forming regions, but 
also along the spiral arms of external galaxies. The sizable sample of extragalactic Cepheids 
detected by the two HST Key Projects aimed at estimating the Hubble constant (Freedman 
et al. 2001; Saha et al. 2001) appears a well-suited extension. The same applies to the large 
sample of Cepheids discovered in M31 and M33 by the DIRECT experiment (Bonanos et al. 
2003, and references therein). 

Finally, we would also like to mention that well-observed open clusters hosting sizable 
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Cepheid samples can be safely adopted to test the accuracy of the physical ingredients 
currently adopted to construct evolutionary models. The difference between the cluster age 
based on the luminosity function and the pulsation ages based on PA and PAC relations can 
supply tight constraints on He-burning phases. Note that these two methods do not depend 
at all upon the cluster distance. 
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Fig. 1. — Theoretical H-R diagram for intermediate- mass stars with different chemical 
positions (see labeled values). See text for more details. 
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Fig. 2. — Predicted Period-Age relations for first-overtone (top) and fundamental (bottom) 
Cepheids with different chemical compositions (see labeled values). The vertical bars display 
the standard deviation of analytical relations. 
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Fig. 3. — Comparison between current Period-Age relation at solar chemical composition 
(solid line) and similar relations available in the literature. The dashed line shows the PA 
relation provided by Magnier et al. (1997), while the dotted and the dashed-dotted lines the 
PA relations by Efremov & Elmegreen (1998) and by Grebel & Brandner (1998), respectively. 
Note that PA relations plotted in this figure have been obtained using different evolutionary 
models. See text for more details. 
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Fig. 4. — Evolutionary ages (Myr) as a function of pulsation period (days) for different 
chemical compositions (see labeled values). Circles refer to fundamental periods estimated 
according to evolutionary predictions that either neglect (Canonical, filled) or account (Non- 
canonical, open) for convective core overshooting. 
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Fig. 5. — Projection onto a plane of Period-Age-Color relations for fundamental Cepheids 
at different chemical compositions (see labeled values). 
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Fig. 6a. — Fit of empirical Color- Magnitude Diagrams in V,B-V (left) and in V,V-I (right) 
of LMC clusters with stellar isochrones at fixed metal content (Z=0.01). Open circles and 
squares mark fundamental and first overtone cluster Cepheids respectively. Individual cluster 
ages, and reddening values are listed in Table 6. 
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Fig. 6b.— 



Same as Fig. 6a but for a different sample of LMC clusters. 
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Fig. 6c— 



Same as Fig. 6a but for a different sample of LMC clusters. 
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Fig. 6d. — Same as Fig. 6a but for a different sample of LMC clusters. 
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Fig. 7a. — Same as Fig. 6a but for a sample of SMC clusters and stellar isochrones with a 
metal content of Z=0.004. 
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Fig. 7b- 



Same as Fig. 7a but for a different sample of SMC clusters. 
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Fig. 7c. — Same as Fig. 7a but for a different sample of SMC clusters. 
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Fig. 8. — Relative difference between evolutionary and pulsation cluster ages as a function 
of isochrone ages for LMC clusters that host at least two Cepheids. The top panel shows 
the difference with the PAC relation, while the bottom one with the PA relation. 
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Fig. 9. — Same as Fig. 8, but for SMC clusters. See text for more details. 
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Fig. 10. — Same as Fig. 6a, but for two Galactic clusters and stellar isochrones with a metal 
content of Z=0.02. 
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Table 1. Input parameters of the pulsation models adopted to compute the PA and the 

PAC relations. 
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a Stellar mass (solar units). ^Logarithmic luminosity (solar units). c Metal (Z) abun- 
dances by mass; the helium abundances are Y=0.25 (Z=0.004, Z=0.008) and Y=0.28 
(Z=0.02). ^Pulsation models constructed by adopting a M/L relation based on evo- 
lutionary models that neglect the convective core overshooting (Bono et al. 1999). 

e Pulsation models constructed by adopting a M/L relation that mimics the behavior of 
evolutionary models that account for the convective core overshooting during H-burning 
(Chiosi et al. 2003; Bono et al. 1999). 
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Table 2. Analytical relations (logT e = a + (3 log L/L & ) for the boundaries of the 
instability strip at different chemical compositions. 
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a Metal abundance. b Zero Point. c Slope. d Standard deviations. 

e Pulsation models constructed by adopting a M/L relation based on evolu- 
tionary models that neglect the convective core overshooting (Bono et al. 1999). 

^Pulsation models constructed by adopting a M/L relation that mimics the behav- 
ior of behavior of models that account for the convective core overshooting during 
H-burning (Chiosi et al. 2003; Bono et al. 1999). 9 Pulsation models from Petroni 
et al. (2003). 
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Table 3. Analytical relations (logT e = a + /Slog L/L Q + 7(logL/L ) 2 ) for the boundaries 
of the canonical fundamental instability strip at different chemical compositions. 
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a Metal abundance. b Zero Point. c Coefficient of 
the linear luminosity term. d Coefficient of the quadratic 
luminosity term. e Standard deviation. f Pulsation 
models from Petroni et al. (2003). 
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Table 4. Period- Age relations (log t a = a + (3 log P) at different chemical compositions. 
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First-overtone 




0.004 


8.4f±0.07 -1.07±0.02 


0.07 


O.Of 


8.29±0.08 -0.80±0.03 


0.08 


0.02 


8.08±0.04 -0.39±0.04 


0.04 




Fundamental 




0.004 


8.49±0.09 -0.79±0.01 


0.09 


O.Of 


8.4f±0.f0 -0.78±0.01 


O.fO 


0.02 


8.3f±0.08 -0.67±0.01 


0.08 



a Cepheid age in years, 
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e Standard deviations. 
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Table 5. Period and ages for selected evolutionary models. 
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3.739 





,870 


97.85 


8 


3.879 


3.716 


1.265 32.36 


4.049 


3.697 


1 


.488 


34.86 


10 


4.163 


3.703 


1.496 19.24 


4.354 


3.678 


1 


.756 


23.56 








Z=0.0198, Y= 


=0.273 










5 


2.967 


3.757 


0.457 100.92 












8 


3.760 


3.693 


1.245 29.93 


3.951 


3.679 


1 


,469 


32.79 


10 


4.073 


3.668 


1.543 18.89 


4.296 


3.648 


1 


.817 


21.63 



a Stellar mass (solar units). ^Luminosity (solar units). c Effective tempera- 
ture (K). ^Period (days). e Stellar age (My). f Metal (Z) and helium (Y) 
abundances by mass adopted to construct evolutionary models. 
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Table 6. Period- Age-Color relations (logt a = a + (3 log P + jCI) at different chemical 

compositions. 





a c 


(3 d 


7 e 






(1) 


(2) 


(3) 


(4) 


(5) 






First-overtone 








0.004 


f8.68±0.06 


-f.f84±0.02 


-2.70±0.33 g 





.06 


0.004 


8.f4±0.06 


-1.17±0.02 


0.65±0.08 h 





.06 


0.004 


8.06±0.06 


-1.16±0.02 


0.64±0.08 i 





,06 


O.Of 


f8.27±0.08 


-0.89±0.04 


-2.63±0.6f g 





.08 


O.Of 


8.04±0.08 


-0.89±0.04 


0.53±0.f4 h 


n 
u 


.uo 


O.Of 


7.94±0.08 


-0.88±0.03 


0.6f±0.f5 i 





.08 


0.02 


14.15±0.04 


-0.48±0.05 


-f.60±0.60 g 





,04 


0.02 


7.94±0.04 


-0.47±0.05 


0.3f ±0.12 n 





.04 


0.02 


7.89±0.04 


-0.46±0.05 


0.36±0.f4 i 





,04 






Fundamental 








0.004 


f4.65±0.09 


-0.88±0.02 


-f.63±0.36 g 





.09 


0.004 


8.34±0.09 


-0.90±0.02 


0.35±0.06 h 





.09 


0.004 


8.24±0.09 


-0.88±0.02 


0.42±0.08 i 





.09 


O.Of 


f6.f3±0.09 


-0.89±0.03 


-2.04±0.42 g 





,09 


O.Of 


8.24±0.09 


-0.90±0.03 


0.36±0.07 h 


0. 


.09 


O.Of 


8.f3±0.09 


-0.89±0.03 


0.49±0.fO i 





,09 


0.02 


ff.60±0.08 


-0.75±0.03 


-0.87±0.34 g 





.08 


0.02 


8.25±0.08 


-0.76±0.03 


0.f6±0.06 h 





.08 


0.02 


8.20±0.08 


-0.75±0.03 


0.2f±0.08 i 


0. 


.08 



a Cepheid age in years. b Metal abundance. c Zero 
Point. ^Coefficient of the period term. Coefficient 
of the color term. Standard deviations. 9 The color 
term is log T e . h The color term is B — V . *The color 
term is V — I. 
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Table 7. Age difference between canonical PA and PAC relations at selected periods. 



logP a 


PA b 


PAC(BE) C 


PAC(RE) d 


PA b 


PAC (BE) c 


PAC(RE) d 


PA b 


PAC (BE) c 


PAC(RE) d 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 






Z=0.004 




Z=0.00l 


I 




Z=0.02 




1 


7.70 


7.64(3.76) 


7.72(3.71) 


7.63 


7.57(3.76) 


7.67(3.71) 


7.64 


7.60(3.73) 


7.65(3.68) 


1.5 


7.30 


7.23(3.74) 


7.33(3.68) 


7.24 


7.21(3.72) 


7.30(3.67) 


7.30 


7.27(3.68) 


7.31(3.64) 



a Selected period. ^Logarithmic age based on the PA relation. c Logarithmic age based on 
the PAC relation along the blue edge of the instability strip. The number in parentheses is the 
logarithmic temperature in which the period is equal to logP = 1, or 1.5. ^Logarithmic age based 
on the PAC relation along the red edge of the instability strip. The number in parentheses is the 
logarithmic temperature in which the period is equal to logP = 1, or 1.5. 
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Table 8. Age estimates for LMC and SMC open clusters that host at least two Cepheids. 



ID a 


b 


E(B-V) C \ogt ISO A 


logtp4 e 


f 

Up A 


log t pAr g 


(j p 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


LMC0394 


2 


0.138 


8.15 


8.015 


0.015 


8.015 


0.005 


LMC0261 


3 


0.142 


7.85 


7.958 


0.033 


7.898 


0.024 


LMC0164 


4 


0.144 


8.08 


7.945 


0.255 


7.870 


0.233 


LMC0164 


3 


0.141 


8.08 


8.069 


0.066 


7.984 


0.054 


LMC0142 


3 


0.152 


7.65 


7.598 


0.156 


7.547 


0.183 


LMC0461 


3 


0.120 


7.85 


7.763 


0.129 


7.795 


0.097 


LMC0457 


7 


0.120 


7.93 


7.737 


0.068 


7.807 


0.073 


LMC0411 


9 


0.145 


7.98 


8.018 


0.075 


7.997 


0.073 


LMC0395 


2 


0.138 


8.08 


8.012 


0.015 


8.015 


0.005 


LMC0622 


6 


0.148 


8.06 


7.996 


0.029 


7.910 


0.019 


LMC0591 


2 


0.147 


8.11 


8.009 


0.006 


7.948 


0.004 


LMC0565 


2 


0.121 


8.06 


8.051 


0.045 


8.019 


0.020 


LMC0559 


3 


0.121 


7.90 


7.946 


0.062 


7.914 


0.046 


LMC0648 


9 


0.175 


7.95 


7.955 


0.079 


7.997 


0.059 


LMC0636 


6 


0.185 


7.85 


7.971 


0.097 


7.887 


0.057 


LMC0633 


13 


0.162 


8.04 


7.967 


0.083 


7.874 


0.078 


SMC0048 


3 


0.101 


8.43 


8.468 


0.184 


8.464 


0.191 


SMC00A8 


2 


0.101 


8.43 


8.368 


0.089 


8.369 


0.138 


SMC0045 


3 


0.094 


8.32 


8.268 


0.036 


8.233 


0.007 


SMC0024 


2 


0.089 


8.28 


8.340 


0.028 


8.291 


0.077 


SMC0016 


5 


0.089 


8.28 


8.284 


0.124 


8.252 


0.142 


SMC0016 


3 


0.089 


8.28 


8.352 


0.067 


8.325 


0.045 


SMC0141 


3 


0.079 


8.22 


8.192 


0.044 


8.138 


0.057 


SMC0105 


2 


0.097 


8.38 


8.368 


0.088 


8.420 


0.185 


SMC0074 


5 


0.094 


8.04 


7.957 


0.222 


7.916 


0.163 


SMC0068 


2 


0.101 


8.30 


8.315 


0.032 


8.270 


0.018 


SMC0216 


3 


0.194 


8.22 


8.199 


0.128 


8.211 


0.177 


SMC0159 


2 


0.084 


8.31 


8.367 


0.119 


8.336 


0.122 


SMC0158 


3 


0.084 


8.00 


8.038 


0.103 


7.972 


0.118 


SMC0149 


4 


0.079 


8.04 


8.111 


0.053 


8.051 


0.064 
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Table 8 — Continued 



ID a 

(1) 


N c b E(B — 
(2) (3) 


V) c log t IS O d \ogtpA° (TPA 

(4) (5) (6) 


log t PAC g 

(7) 


<JPAC 

(8) 





a Cluster identification according to Pietrzynski & Udalski (1999). 
b Number of cluster Cepheids. 

c Mean cluster reddening based on individual cluster Cepheid red- 
denings (Pietrzynski & Udalski 1999). 

d Cluster age based on stellar isochrones (yr). The uncertainty 
on this estimate is of the order of a IS o = 0.075 dex. 

e Mean cluster age based on fundamental and first overtone PA 
relations (yr). 

f Standard deviation. 

g Mean cluster age based on fundamental and first overtone PAC 
relations (yr). 

h Standard deviation. 
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Table 9. Age estimates for two Galactic open clusters that host at least two Cepheids. 



ID a 

(1) 


logP b 

(2) 


< B > c 
(3) 


< V > c 
(4) 


</> c E(B-V) d \ogt PA c logfgic V)f 
(5) (6) (7) (8) 


iOg IpAC 

(9) 










NGC 7790 h 








CEa Cas 


0.711 


12.070 


10.920 


0.562 


7.834 


7.788 




CEb Cas 


0.651 


12.220 


11.050 


0.548 


7.874 


7.827 




CF Cas 


0.688 


12.335 


11.136 


9.754 0.531 


7.849 


7.824 


7.838 










NGC 6067 h 








QZ Nor 


0.578 


9.774 


8.866 


7.893 0.249 


7.923 


7.872 


7.865 


V340 Nor 


1.053 


9.526 


8.375 


7.151 0.315 


7.605 


7.704 


7.708 



a Cluster variable identification. ^Logarithmic Period [days]. c Mean B,V,I magnitudes 
(mag). ^Reddening. e Age estimate based on the PA relation (yr). f Age estimate based 
on the PAC(b_vo relation. 9 Age estimate based on the PAC(y_/) relation. ^Photometry 
and reddening estimates for cluster variables from Berdnikov (2000) [NGC 7790], and Laney 
& Stobie (1994) [NGC 6067]. 



